clear;clc;
A=[9,5;2,6];
x=randi([1,500],2,1);
[P,D]=eig(A)
syms n
syms x_n
P*D^n*P^-1  %求出A^n
x_n=ans*x  %求出xn
limit(x_n(2,1)/x_n(1,1),n,inf)
x